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Transport of massless Dirac fermions in graphene monolayer in the presence of a combination 
of singular magnetic barriers and applied electrostatic potential is analyzed. The similarity of 
such transport with the transmission of light through a medium with modulated refractive index 
discussed for singular magnetic barriers recently^ has been extended for the case of an added 
electrostatic potential. We have also calculated the quantum version of the Goos Hanchen shift that 
the electron wave suffers upon being totally reflected by such barriers. Additionally, we have found 
that the presence of such electric and magnetic barriers modifies the band structure of graphene 
substantially near the Dirac point. This affects the transport near the Dirac point significantly and 
has important consequences for graphene based electronics. 
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I. INTRODUCTION 

Charge carriers in monolayer graphene behave as mass- 
less Dirac fermions leading to an intriguing set of trans- 
port phcnomenaii 2 ^^^. One of the consequences of 
this Dirac fermion-likc behaviour is Klein tunneling 7 -^, 
where an electron can pass through a high barrier in 
contrast to the conventional tunneling of non-rclativistic 
electrons. This has consequences which are exper- 
imentally observable^. Recent experiments demon- 
strated Klein tunneling in graphen o 11 ' 12 and in carbon 
nanotube*^. Crucial to the designing of graphene based 
electronics is to attain confinement of electrons within 
a mesoscopic or nanoscopic size of the sample. Due to 
Klein tunneling, it is difficult to confine electrons conven- 
tionally. For this reason, several alternatives have been 
suggested. In a recent work it has been suggested that 
a magnetic barrier can effectively block Klein tunneling 
and achieve confinement for such massless Dirac fermions 
for graphene^. This has motivated a number of subsc- 
quent reports^^^^^^* 2 ^ 2 ^. 

In the ballistic regime, the scattering of electrons 
by potential barriers can be understood in terms of 
phenomena like reflection, refraction and transmission, 
leading to an analogy between electron transport and 
light propagatio n 26 ' 27 . This provides better insight in 
both fields and often leads to novel device applica- 
tions. For example, the similarity between transmis- 
sion of de Broglie waves of two dimensional electrons 
satisfying the Schrodinger equation in the presence of 
a one-dimensional electrostatic potential barrier and the 
propagation of light is well established and can be used 
for lensing and focusing of electron o 28 ' 29 . Light prop- 
agation through optical fibers can also be understood 
in terms of the Schrodinger equation through various 
potential barriers 3 ^. Extension of such optical analo- 
gies to understand transport in graphene is character- 
ized by the massless Dirac fcrmion-like nature of charge 



carriers in graphene. This was clearly demonstrated by 
Cheianov et a& by showing that electron transport in 
graphene with an applied split gate voltage is akin to 
light propagating through a metamaterial with negative 
refraction inde x 32 ' 33 . In a subsequent work^i, it has 
been shown that such electron propagation can be in- 
terpreted in terms of the quantum version of the Goos- 
Hanchcn effect. The possibility of guided modes in a 
graphene waveguide was proposed recently using similar 
arguments 3 -*. Guided modes through graphene constric- 
tions in the sub-wavelength regime of coherent transport 
using optical analogy was also studied recently^ 6 -. 

Most such analogies with optical phenomena have been 
drawn only for electrostatic potentials. Given the im- 
portant role magnetic barriers play in confining mass- 
less Dirac fermions, a similar analogy is very much 
needed for the case of electron transport through mag- 
netic field. The approach used to establish an optical 
analogy for electrostatic potentials cannot be easily car- 
ried forward to the case of electron transport through a 
magnetic field. The magnetic field bends the trajectory of 
the electron continuously which exhibits cyclotron mo- 
tion. Thus, a direct analogy with the linear propaga- 
tion of light is not possible. This consideration, however, 
changes when one considers the motion of an electron in 
the presence of a highly inhomogencous magnetic field. 
Particularly if the range of inhomogeneity is much smaller 
than the typical radius of the cyclotron, one is left with 
plane wave like scattering states. A limit of such highly 
inhomogeneous magnetic field is the singular magnetic 
barrier, which gives rise to step function like magnetic 
vector potential (MVP) barriers. Following the work by 
F. M. Peeters et a/.—, this has been widely studied in 
the literatur e 38 ' 39 which shows that such a singular mag- 
netic barrier can lead to wavevector filtering. In a recent 
paper—, it has been shown by us that, for massless Dirac 
fermions in single graphene layer propagating through 
a series of such singular delta function type magnetic 
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barriers^!, one can build an optical analogy and fcrmions 
behave like a light ray passing through an optical medium 
with a modulated refractive index-2ii. However, the corre- 
sponding Snell's law is very different from that in classical 
optics and is not specular. In fact, the analogy can be 
easily extended to the case of ordinary non-relativistic 
electrons. Changing the direction and magnitude of the 
magnetic field of such barriers changes the refractive in- 
dex, and this generates a number of interesting phenom- 
ena that have direct analogues in optics. For example, it 
is possible to explain the absence of Klein tunneling and 
the resulting confinement^ in presence of such magnetic 
barriers using the language of optics. This analogy can be 
extended to propose practical structures like Bragg reflec- 
tors and resonant cavities using such magnetic barriers. 
Similar analogy has also been recently studied by com- 
paring the electron transport in periodic or disordered 
array of smooth scalar potentials in graphene supcrlat- 
tices with the propagation of electro-magnetic waves in 
a layered medium^!. 

In this paper, we explore the transport of massless 
Dirac fcrmions in graphene under the combined effect 
of a magnetic barrier and an electrostatic voltage such 
as the one used in Ref.— to explore an exotic transport 
regime. The purpose is two fold. One, in graphene, 
electrons behave as massless Dirac fermions and have 
a linear band structure, albeit only close to the fermi 
level. In this region, small electrostatic potentials can 
greatly affect the electronic states by shifting the fermi 
level and causing the graphene sheet to behave as cither 
an electron-deficit (p-type) or a hole-deficit (n-type) ma- 
terial. For this reason, it is important to include the ef- 
fect of electrostatic potentials on any proposed structures 
made with graphene. The main findings in the paper can 
be summarized as follows. Firstly, the electron trans- 
port through a number of combinations of electrostatic 
and magnetic barrier is explained using the language of 
geometrical optics. This discussion brings out certain 
characteristics of transmission like high asymmetry, the 
non-specular nature of the refraction of electron waves as 
well as high rectification properties. As a byproduct of 
this optical analogy, we describe the Goos-Hanchen shift 
of massless Dirac fermions at such barriers which can 
change sign as well as magnitude abruptly upon total in- 
ternal reflection at various incident angles. An analysis 
of transport through a finite and an infinite series of such 
EMVP barriers concludes our investigation. This treat- 
ment also provides a way to control the band structure 
as well quasiparticle transport in monolayer graphene. 



II. ELECTRON PROPAGATION AND ITS 
SIMILARITY WITH OPTICS 

To begin with, we shall briefly review the optical 
analogues of non-relativistic and rclativistic electron 
transport through various electrostatic as well as mag- 
netic barriers. Consider a two-dimensional electron gas 



(2DEG) that is now routinely produced using semicon- 
ductor heterostructures. When an electron at the fermi 
level Ep is incident on a potential barrier V, its momen- 
tum parallel to the interface of two regions 1 and 2 is 
conserved; i.e., pisin^i = p2sin#2, where pi t 2 are the 
momenta and are the angles in the two regions. This 
leads to the following Snell's law for electrons^ 

^r = a-^ (D 

sm#2 Ef 

When a similar scattering problem is considered for 
the two-dimensional case of massless Dirac fcrmions in 
graphene, in the absence of any potential barrier the 
fermi level lies at the Dirac point where the conduction 
and valence band touch each other. Thus the electronic 
states near the fermi level are equally populated by elec- 
tron and holes and the system is charge neutral. The 
introduction of potential barrier locally raises the fermi 
level from this charge neutral Dirac point by placing 
it inside the conduction band. Thus inside the barrier 
regime electron states near the fermi level are now pop- 
ulated by electrons only and it becomes electron reach 
(n-type). However, because of the chiral nature of such 
massless Dirac fcrmions, this conversion of charge carriers 
from hole to electron is accompanied by a corresponding 
change in the direction of the momentum vector. Thus, 
the equality of the momentum component parallel to the 
interface gives rise to a negative refractive index^I. 

If a similar Snell's law for electron transport in the 
presence of a magnetic field is to be achieved, we need 
to find a suitable magnetic field profile which will scatter 
the electrons in much the same way as an electrostatic 
potential docs. As discussed earlier, this implies that the 
electron cannot be allowed to complete the cyclotron tra- 
jectory. Consequently, the magnetic field has to be highly 
inhomogeneous on the scale of the fermi wavelength. 

An extreme case of such an inhomogeneous magnetic 
field is the one introduced in Ref.— having the following 
profile of the transverse magnetic field B and the corre- 
sponding vector potential A in the Landau gauge: 

B = B z {x)z = B£ B [S(x + d) - S(x - d)]z 
A v {x) = B£ B e{d 2 -x 2 )y (2) 

Here £b = VhceB is the magnetic length. Such a mag- 
netic field creates a wavevector dependent potential bar- 
rier which scatters electrons, and such a magnetic vector 
potential (MVP) barrier can be used for wavevector fil- 
tering. It is thus a very natural question to ask whether 
electron scattering by such MVP barriers admits an op- 
tical analogy very similar to the one suggested through 
Eq. ([1]). The problem attains additional importance in the 
case of massless Dirac fermions in graphene since Klein 
tunneling does not occur through such inhomogeneous 
magnetic barriers leading to confinement of electrons^. 
This problem has been investigated recently^ for the 
case of massless Dirac fcrmions in graphene and it has 
been shown an optical analogy can indeed be built. But 
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the resulting Sncll's law is exotic since it is non-specular 
in nature unlike similar laws in optics as well as for the 
motion of an electron through electrostatic potential bar- 
riers. Namely, the matching of the momentum compo- 
nents and the energy conservation when applied to the 
problem of electron scattering in the presence of magnetic 
barrier like the one given in Eq.@ gives 
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Here, [k x , k y ] = kp[cos<f>, sin<f>] outside the magnetic bar- 
rier and 4> is the incident angle for an electron wave and 
[Qx> k y — j^] = fcp[cos 0, sin 0], where 9 is the angle of re- 
fraction. The relation given in Eq.Q implies that, for a 
wave incident with positive <j), the wavevector will bend 
towards the normal. Similarly, for a wave incident with 
negative incidence angle, the corresponding wavevector 
will bend away from the surface normal inside the bar- 
rier region. A scries of such magnetic barriers will thus 
lead to highly asymmetric transmission of electron a 16 ' 17 . 
Eq.(|3J) also yields that when | sin \9\\ > 1 the angle 9 be- 
comes imaginary and total internal reflection(TIR) oc- 
curs. This naturally explains why such barriers can con- 
fine electrons. Confinement will occur when sin|0| > 1 
for -f < 4> < and when sin \9\ < -1 for < 4> < f ■ In 
the latter case, this requires the wavevector to be nega- 
tively refracted at sufficiently high magnetic field before 
TIR occurs. 

It has been pointed out^ 7 - that, to obtain symmetric 
transmission, one needs to place two such single MVP 
barriers side by side but orient them oppositely. The 
resultant magnetic field in Landau gauge is given by 
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B z {x)z = B£ B [S(x + d)- 2S(x) + S(x - d)]z(4) 



The energy conservation in medium 1 (— d < x < 0) and 
medium 2 (0 < x < d) leads to 
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The incident angle is — \ < <f> < \ and the angles of 
refraction are 0i and 9 2 in media 1 and 2 respectively. 
Thus, for such double MVP barriers, the wavevector 
bending towards (away from) the surface normal in the 
first half of the barrier bends away from (towards) the 
surface normal in the second half of the barrier achiev- 
ing symmetric transmission. Because of total internal 
reflection from the first as well as second half of a double 
MVP barrier beyond a critical angle of incidence, the re- 
flectivity of such a barrier is relatively higher than that 
of a single MVP barrier. As shall be discussed later on in 
this paper, a series of such MVP barriers put side by side 
can work as a Bragg reflector, and the associated band 
structure also shows the effect of the magnetic field on 
the transport properties. 



III. TRANSPORT IN A COMBINED MVP 
BARRIER AND ELECTROSTATIC POTENTIAL 

The electron focusing property of a single layer of 
graphene in the presence of a split gate voltage due to 
negative refractior^i has already been identified. In this 
section, we shall see how this behaviour gets modified in 
the presence of additional MVP barriers. We shall con- 
sider how electrons will behave when one has a potential 
step (created by a local gate voltage) on top of an MVP 
barrier. We shall discuss this for a number of combina- 
tions of electrostatic and MVP (EMVP) barriers which 
we consider simple but illustrative. 



A. Single MVP barrier with electric field 

Let us begin by considering the transmission of mass- 
less Dirac electrons through the following combination of 
a single MVP barrier described in Eq. @ and a commen- 
surate electrostatic step potential barrier, namely 



V(x) = V,A(x) = B£ B , 
V(x) = 0,A(x) = 0, 



|a;| < d 
\x\ > d. 



(6) 



We assume V > and discuss the case V < later on. 
For definitcness, we shall set the energy of the incident 
electrons as E — Ep, where Ep is the fermi level located 
at the charge neutral Dirac point in single layer graphene. 
The voltage V is measured with respect to this Ep at 
the Dirac point. However, the entire treatment that will 
follow is applicable in the neighbourhood of charge neu- 
tral Dirac points as long as the dispersion remains linear. 
Additionally, we also assume that scattering by such an 
EMVP barrier is not so strong that it will break the de- 
generacy of the K and K' points. The typical length 
scales over which such a delta function approximation of 
the magnetic barrier is valid has been detailed in ReL^ 7 -. 

Following discussion in section HT1 the electrostatic po- 
tential barrier shifts the fermi level from the charge neu- 
tral Dirac point either in conduction band for V > or 
in the valence band for V < 0. This creates locally an n— 
or p— doped region. Thus, the following analysis can be 
used to understand the transmission through such n or p 
regions in the presence of an MVP barrier. We shall call 
these regions electrostatic plus magnetic vector potential 
(EMVP) barriers. In the region — d < x < d at either of 
the K points, the motion will be described by 



Vp 



v 

Tl X + Illy 



7T X — 27T. 

v_ 

VF 



E 



tpi 

>P2 



(7) 



Here 7r = p + and V implies qV where q = — |e| 
Since the vector potential is in Landau gauge, the sta- 
tionary solutions can be written as 

1>i,2(x,y) = (t>i, 2 {x)e ik » v 



4 



V 7 >V 6 >V 5 >V 4 >V 3 >V 2 >V 1 >0 



.•• v=o 




where V = is the potential barrier in the dimcn- 
sionless form. Let us define the angle of incidence of the 
electron wave outside the barrier region as 6 = tan _1 (f^) 

and the angle of refraction 9 = tan -1 ( k q g^ j B _ 1 ) which the 
wavevector made with the surface normal in the barrier 
region. In terms of these angles, Eq. (fT2|) can be rewrit- 
ten to obtain the Snell's law analogue for electron waves 
of such massless Dirac fermions incident on the EMVP 
barrier as 



V 3 >V 2 



FIG. 1: Ray diagram of electron propagation through a sin- 
gle EMVP barrier. With increasing applied voltage, the re- 
fracted angle 02 increases continuously in the clockwise direc- 
tion. Dashed rays are negatively refracted waves. 



Substituting these solutions in Eq.®, one gets the fol- 
lowing coupled one-dimensional equations 
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The above equations can be decoupled to yield 
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The corresponding stationary solutions (j>i t 2{x) in the 
barrier region are similar in form to those for the cases 
of a pure magnetic barrier— or an electrostatic step 
potential^ and are given by 



e ikxx _|_ rp ik x x jf x < _^ 
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te lkxX if x > d 
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Here {k x , k y } and {q x , k y — j^} respectively refer to the 
x and y components of the wavevector inside and outside 
the barrier regime. Since the incident energy E — Ef = 
Hvpkp, substitution of solutions (fTT)|) and (fTTj) in Eq.© 
leads to 
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Comparison with the expression in Eq. Q shows that the 
effect of the potential barrier is therefore to scale the 
refraction angle by the factor E ^_ v . In the present case, 
E = Ep- Thus, for < V < Ep, this scale factor varies 
from 1 to oo. For V > Ep, the scale factor becomes 
negative and its value increases from — oo to as V is 
increased to oo. We shall now discuss how the incident 
wave will get refracted in the barrier regime due to this 
non-monotonic and discontinuous behaviour of the scale 
factor. 

For definiteness, let us first consider the cases when 
< (f> < Also, we use not too high a magnetic barrier 
field strength (B) such that sin#|y = o = sin0— fc ^ > 0. 
For V < Ep, the scale factor is greater than 1. This gives 

sin(9(V" > 0) > sin6»|y =0 ,0 < V < E F (14) 

The above relation implies that the angle of refraction 
will start increasing with V . As a result, the wavevector 
will bend increasingly away from the surface normal. For 
V = 0, sin 9 < sin<p , suggesting the electrons are going 
from a rarer to a denser medium. With increasing V, 9 
will increase for a constant <f> and electrons will behave 
like passing into an increasingly rarer medium. At some 
point, sin# will be greater than sin</>, making the barrier 
regime behave like a rarer medium as compared to the 
region outside the barrier. If V is such that 
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the electron wave will suffer total internal reflection at 
this junction. Thus, the following result implies that for 
a given strength B of MVP barrier and a given angle of 
incidence <f> > sin" 1 fc ^ , by increasing the height of the 
electrostatic barrier, it is possible to totally reflect the 
electron wave from the barrier. 

The point V = Ep represents a singularity in the spec- 
trum and demands a separate discussion. In the ab- 
sence of magnetic barriers, such a point represents the 
zero modes for Dirac operators and leads to the emer- 
gence of new Dirac points. This has been discussed in a 
number of recent works considering Andreev Reflection 
in a graphene based NIS junction^, for the electronic 
states of monolayer graphene monolayer in a periodic 
potcntia l 43 i 44 ' 45 . In presence of a magnetic barrier these 
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zero modes become the zero modes of the of the modi- 
fied Dirac operator. The nature of corresponding solu- 
tions are different. The magnetic barrier breaks the time 
reversal symmetry explicitly. The equations satisfied by 
<f>i 2 will be 



B=0.1T 



-dl 



fl,2 







Thus, the solutions along the x-direction are either ex- 
ponentially decaying or growing. Since these solutions 
exist in the region of total internal reflection we shall re- 
tain only the decaying one. This can be contrasted with 
the case without any magnetic barrier— where the equa- 
tion obeyed by the zero mode solution is —d%4>i,2 = 0. 
yielding linear solutions instead of an exponential one. 
The significance of such zero modes is that on the twe 
sides of the singular point V = Ep, the relative sign be- 
tween the angle of incidence and the angle of refraction 
becomes opposite. 

Beyond the singular point at V = Ep, the scale factor 
becomes negative and thus the sign of the refraction angle 
becomes opposite to the sign of the angle of incidence, 
implying negative refraction. In this case, the region with 
EMVP barrier will act like a left-handed metamaterial. 
In this negative refraction region, as long as 
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the electron wave will remain total internally reflected. 
However, at very high values of V, 
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'Ep-V 
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(15) 



The electron wavevector will again retrace its path back 
to the first medium by negatively refracted. As V — > 
oo, sin 9 — > 0. If, on the other hand, > <j) > — f i 
for V = 0, sin \9\ > sin |0| . Thus, electrons incident 
with such an angle on a pure MVP barrier are seen as 
passing from a denser to a rarer medium. Now, if the 
electrostatic voltage is turned on and < V < Ep, the 
angle of refraction \9\ will increase with V and eventually 
the electron wave gets totally reflected internally. Thus, 
again by increasing V it is possible to totally reflect an 
electron wave for any given (f> and B. For V > Ep, 
again the sign of sin 16*1 will be opposite to the sign of 
sin|(/)|, thus making the refraction negative. For very 
high value of V, sin |6*| will again retrace its path back in 
the barrier regime, but since its sign is opposite to the 
sin|f3| the electron wave gets negatively refracted. This 
situation is depicted in Fig[T] where we have shown how 
the wavevector for the refracted ray changes with the 
increasing value of the voltage V. 

The preceding discussion is summarized in FigJ5] by 
plotting the refraction angle as a function of the angle of 
incidence and the strength of the electrostatic potential 
barrier respectively for B — 0.1 Tcsla and B = 3 Tesla. 
Along the y or axis, V = and the scale factor is 1. 



denser->rarer 
rarer->denser 



0.5 






gap 


negative refraction 


rarer->denser 
denser->rarer 



1.5 

qV/E p 



FIG. 2: Phase diagram of the refraction angle at magnetic 
field 0.1 Tesla and 3 Tesla. The gap corresponds to the region 
where total internal reflection (TIR) occurs, d — lOOnm for 
this as well as for all subsequent figures. 



This is the case of MVP barrier without any electrostatic 
potential. Since (f> goes from to sin 9 goes from 

"d^ t0 1 " W For °- 1T ficld in Fi sH a )> Eq.® 
yields a real solution for 9. Now, for a given value of <ft 
such that sm9 is positive on the c/)-axis, along the x or V- 
axis the scale factor changes in the way described in the 
preceding paragraph. Accordingly, sin# passes through 
the following four phases: (a) sin# < sin cf> < 1 (rarer 
to denser); (b) 1 > sin 6* < sine/) (denser to rarer); (c) 
| sin 9\ > 1 (TIR (gap)); and (d) < sin6» < -1 (negative 
refraction). Thus, in the (c) region, 9 cannot be real 
and that implies total internal reflection. In this region 
where total internal reflection occurs, either sin# > 1 or 
sin 9 < — 1 and these two regions are separated by the line 
Ep = V . In the negative refraction region, the system 
will behave like a metamaterial. 

If the incidence angle <fi is smaller than s\oT 1 { k ^ gB ), 
sin 6? is negative on the y axis. The corresponding sit- 
uation has been depicted in the lower part of the plot. 
Since the scale factor along the V axis changes in the 
same way, here with increasing V sin 9 will pass through 
(a) a negatively refracted region, (b) a TIR gap region, 
(c) a denser to rarer interface, and (d) a rarer to denser 
interface. Thus, the behaviour in this region changes in 
a reverse way as compared to the preceding case with the 
increasing V . All these different regions meet at a limit- 
ing point given by sin <j> — k ^ and V = Ep where the 
behavior is singular and is described by the zero modes 
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of the Dirac operator in presence of the EM VP barriers. 

All the above mentioned properties will also affect the 
transmission through such a barrier. The transmission 
amplitude can be obtained from Eqs. (fit))) and (fTTj) by 
using the continuity of the wave functions at x = —d and 
x = d as 
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where D = 2d and now q x = kj cos 6. The expression 
of the transmission is same as the one for electrostatic 
potential barrier— or single MVP barrier—. 

The transmittance (T = t*-t) and reflectance R = 1—T 
are of the same form as given in Refi^. In Fig[3l for two 
given strengths B, the transmittance as a function of 
angle of incidence 4> is plotted for a range of values of 
V. The first row plots the transmission coefficient for a 
smaller value of B = 0.1 Tesla. The central figure cor- 
responds to the case for pure magnetic barriers, namely 
when V = 0. We have already noted that for V = Ep, 
since the scale factor E ^ v diverges, the barrier becomes 
fully reflecting. And, the nature of the solutions is differ- 
ent here and corresponds to zero modes. Therefore it will 
be interesting to study the transmission properties as this 
singular point is approached. To see this, V is increased 
to qV = 0.84EF (immediate right from the central plot). 
As one can see, in this case transmission takes place over 
a very small window along the (f> axis asymmetrically lo- 
cated in one quadrant of the the <fi axis. In the immediate 
lower plot where B is changed from 0.1T to 3T, a sim- 
ilar narrow window of transmission gets more shifted to 
one side of the d> axis. This is because the deviation 



from the normal direction 



increases with B. On 



the other hand, if V is changed to the opposite sign such 
that qV = — 0.84Ep, the scale factor becomes E ^ v \ , 
which is always < 1. Thus, the angle of refraction on 
both sides of the surface normal will now be decreased 
upon multiplication with this scale factor as compared 
to the case when V = 0. As a result, the critical angle 
for total internal reflection will go up. In the plot im- 
mediately left of the central column, we see that for the 
same magnitude but of opposite sign of V transmission 
exists over a large range of <fi. If \V\ is increased to 4.2Ep 
such that |V| 3> Ep, the barrier becomes transmitting 
for all values of <fi. This can be seen from the extreme 
right and left plots. This also agrees with the features of 
the refraction map in Figf2^a), where at higher values of 
V along most of the region of the </> axis the refraction 
angle 9 is real. 



in the transmission in presence of a pure magnetic bar- 
rier (when V = 0) increases with B. To summarize, the 
above analysis clearly establishes that, for not so high V, 
the transmission through the barrier can be drastically 
changed by switching V from positive to negative. Thus, 
such a barrier at certain values of V possesses important 
rectification properties. 



B. Double MVP barrier with split gate voltage 

In this section, we shall consider transmission through 
a double MVP barrier added with a split gate voltage, 
henceforth called a DEMVP barrier. The vector and 
scalar potentials are characterized by 

V(x) = -V, A y (x) = B£ B , -d < x < : region I 
V(x) = V 7 A y (x) = ~Be B ,Q<x<d: region 11(17) 

In regions I and II, the local fermi level Ep is respec- 
tively lowered and raised from the charge neutral dirac 
point to valence band and conduction band. Thus the 
region / and II becomes p and n type region. Thus, the 
split gate voltage creates an p — n junction on monolayer 
graphene. A reversal in the sign of the voltage can simi- 
larly create an n — p junction. This has interesting trans- 
port properties that can be used for the purpose of elec- 
tron focusing^. In the following analysis, we shall learn 
about the effect of highly localized magnetic barriers on 
the transport through such ann-p junction. Assuming 
that the fermi level lies very close to the charge-neutral 
Dirac point, the transport electrons obey 



±- 



■ iif„ 



±- 



= E 



02 



(18) 



where the — ve sign is for region I and the +ve sign is for 
region II. For definiteness, we again set E = Ep albeit 
the treatment is valid in the entire linear spectrum region. 
The stationary solutions in both these regions in Landau 
gauge are again of the form ^1,2(2;) = <t>i,2{x)e lky , where 
the subscripts 1, 2 denote regions I, II respectively. Sub- 
stitution of these solutions in the equations l|18p respec- 
tively gives in region I (— d < x < 0) the equations 



d 1 

hv F [-i— ~i(k y - — )]</> 2 = {Ep + V)^! 
ox Ib 

9 1 
Hv F [-i— + i(k y - — )}(/>i = (E + V)fo (19) 
ox Ib 



In the lower row of FigJ3] is a similar plot but for a 
higher value of B = 3T. Since the higher magnetic field 
also means a much higher reflectivity of the barrier, the 
effect of turning on an electric field is to introduce more 
asymmetry in the transmission as one compares T(cf>) 
with T(— <fi). This is due to the fact that the asymmetry 



Similarly, the equations in region II (0 < x < d) are 

1 



d 

hv F [-i— - i(k v + —)}<t> 2 = {Ep-V)^! 



KM-i-^ + + j-)]<h = {E-V)4> 2 (20) 
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FIG. 3: Transmittance T through a single EMVP barrier at B = 0.1 Tesla and 3 Tesla at different angles of incidence cj>. 



V 





Bz(x) 







-d 




B,(a?) 


d 




-V 







FIG. 4: Electric and magnetic potentials for a DEMVP bar- 
rier. 



Each pair of equations can now be decoupled to yield 
d 2 „ l_ 



(hv F ) 2 [ 
(hv F ) 2 [ 



dx 2 
dx 2 



■) 2 ]0i,2 = {Ef + V) 2 ^, 



(ky 



Ib' 



if - d < x < 
ii,2 = (E F -V) 2 <t> lt2 , 



The stationary solutions that satisfy the above equations 
are left and right moving plane waves of the form 



»1,2 OC 



e lqix , if - d < x < 

e iq 2 x if < x < d 



(22) 



Substituting in the above equations and setting Ep = 
hvpkp, we get 



*2 



k F , if \x\ > d 



Q 2 i + IK - 

*-B 

<& + [ky + H 



(23) 

if - d < x < (24) 
^} 2 , ifO<x<d (25) 



The magnitudes of the modified wavevectors in regions I 
and II are denoted by hp 2 = ^§^r^ and they are related 
to the incident wavevector kp by 



kU B 



kl£ B 



k F e B + v 

k F i B - V 



(26) 
(27) 



Here, V = j^—- The explicit solutions in the various 
regions can then be written as 

Jtk- 



e x + re 

ae i( i lX + be 
ce iQ2X , de 



■ik x x 
-iqix 
-iq 2 x 



if x < —d 
if - d < x < 
ifO < x < d 



(28) 



s\[ae 
S2\ce 



s [ e i(k m x+<l>) _ re -i(fc^+0)] }f x < _ d 

i(q 1 x+0 1 ) _ he -i{q 1 x+8 1 )^ if < a; < 

i(q2X+0 2 ) _ de -i{q 2 x+e 2 )} ifO < x < 

ste i(k*X+<t>) [f x>d 



tfO<x<d (21) Here, tan 6»i 



(29) 

and tan 62 = — — — , and the angle 

qix qix 

ci, 2 gives the angle between the propagation vector and 
x-axis in medium I, II. Also = sgn(E F ± V) and can 
be ±1 depending on the value of V for a given E F . Using 
the definition of 6>i, 2 and Eqs.|23|) - ([27]) 



sin t>i 



sin 6^2 



g[sin0-sgn(^)^i 



r sin — sgn(0) 



E F + V 
k F 



k F li 



4 [sin0 + sgn(0) fcj 



Ft B 



E F . , 1 
[sin + sgn(0) 



(30) 



(31) 
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Thus, the expression for sinfli^ for the case of a double 
MVP barrier (V = 0) gets multiplied by the scale factor 
E E ± V to yield the expression for sin #^2 for a double 
EMVP barrier, namely when V 7^ . The two scale 
factors behave in different ways. Let us first consider the 
case V > 0. For E F > V, 1 < sS^tt < 00 and for 

V > E F , -00 < s^-r? < 0. 



however, satisfies 1 > 



E F -V 

The other scale factor, 



E+V 



> 0. 




qV/E F 




-0.5 



qV/E c 



FIG. 5: color onhneThe refraction map for a double EMVP 
barrier of field strength (a) 0.1 Tesla and (b) 3 Tesla. The 
explanations for the signs and the colors are given in the text. 

When V < this behaviour gets reversed. All these 
effects have been plotted as a refraction map in Fig|5] In 
this figure is plotted side by side sin 0^2 given by relation 
([50)) and (pTIj) over a range of values of V positive as well 
as negative. The upper and lower figures correspond to 
B = 0.1T and 3T respectively. The white region corre- 
sponds to where one of sin 0^2 becomes imaginary and 
as a result the electron wavevector gets totally internally 
reflected cither from region I or from region II. Such TIR 
occurs when either | sin 9\ | or | sin 62 | becomes > 1 . The 
figure shows that, for a given <fi and B , both | sin 6*1,2 1 
cannot be simultaneously > 1, a condition which follows 
from Eqns. (P01) and (j5T|) . 

The colored portion belongs to where #1,2 are real, and 
for the corresponding values V and 4> there will be trans- 
mission. On the colored portions we have marked the 
sign of the angle of refraction in region I and region II 
in an ordered way. The first sign corresponds to the 



sign of sin \9i \ and the second sign corresponds to sin \02\- 
As can be seen, the transmission probability greatly in- 
creases with the increased value of |V|. This is because 
the relative effect of the electrostatic potential barrier as 
compared to magnetic barrier increases in these regimes. 
Particularly, the plots clearly show that over a range of <fi 
the angles in one of regions I or II which was imaginary 
for a given V becomes real upon increasing \V\. However, 
its sign becomes opposite to that of (f> implying negative 
refraction. This phenomenon explains Klein tunneling 
for high V using the language of optics^. The various 
sectors of the refraction map also reflect the symmetry 
properties of such barriers as one makes the transforma- 
tion V -> -V and <j) — > -<j) in Eqns. (J3H1) and (|3"Tj) . 

One of the consequences of the different behaviour 
of the scale factor in Eq. (|50|) and in Eq.([2I|) is that 
the transport through such double EMVP barriers is 
asymmetric as compared to the symmetric transmission 
through the double MVP barrieril. This is shown in 
FigJS] by plotting the transmission coefficient as a func- 
tion of <p. The transmission can be calculated by gener- 
alizing a transfer matrix approach to calculate the trans- 
mittance through double MVP barriers. We introduce 



the notation A = e 
following matrices: 



-ik x d 



Bi, 



e -»?i*,2*d to define the 



M A 



A 

A* 



1 

S,Sl,2 



,M Bl 



B U 2 
B 



1.2 



(32) 
(33) 



The continuity of the wave functions at x = —d,0,d re- 
spectively yield 



M S M^M A 
M sl M 6l 



M S2 M d2 M% 



M Sl M 6l M Bl 



M S2 Me 2 



M s M<j,M A 



(34) 



The above set of equations can be combined to elimi- 
nate a,b,c,d and to obtain transmission coefficient t and 
reflection coefficient r as 



MXiMsM^TEMVpMjM+MZ 



(35) 



Here, Temvp is the transfer matrix through such double 
barrier and is given by 



T, 



EMVP 



M Sl Mg 1 M Bl (M Sl Mg 1 ) ~ 1 M s . 
M 02 M B2 (M S2 M S2 )- X 



(36) 



The transmittance given by T = t-t* is plotted in Figj6] 
The central column of Fig JS] corresponds to V — 0, where 
we plot transmission through the pure double MVP bar- 
rier with symmetric transmission. At small values of V 
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0.1T, qV/E^-4.2 0.1T, qV/E^-1 .4 0.1T, qV/E^O 0.1T, qV/E^I .4 0.1T, qV/^4.2 




FIG. 6: Transmittance T through a double EMVP barrier for 0.1(upper row) and 3 Tesla (lower row) at different angles of 
incidence <j>. 



the transmission gets asymmetric. In fact, this asymme- 
try is higher for larger B (the lower row of the figure). 
This can be explained by noting that for larger B the 
shift in the longitudinal component of the momentum on 
two sides is more asymmetric. In case of the pure double 
MVP barrier-! 7 ., since the refraction properties of right 
and left side of the barrier get interchanged, they finally 
get compensated and the resultant transmission through 
such a barrier becomes symmetric. However, here in re- 
gion I and region II they get multiplied by different scale 
factors, namely E E ± V ■ Thus, the difference between the 
shift in the right and left moving wavevectors becomes 
larger instead of getting compensated while getting into 
region II. This difference further increases with increasing 
value of the magnetic field. At very high value of the volt- 
age V , however, the transmission characteristic is domi- 
nated by the effect of electric field. Finite transmission 
takes place at all values of incident angle — £ < <j> < £ , 
and this can also be seen from the corresponding refrac- 
tion map in Fig|5l However, the amount of transmission 
changes substantially over the entire range of <f>. As a re- 
sult, at higher values of \V\, the transmission properties 
get less asymmetric (the plots in the extreme right and 
left column). 

This transfer matrix formalism can now be extended 
easily through a sequence of n such barriers which can 
be used to create a Bragg reflector. We have already 
analyzed such structures in our earlier work with pure 
magnetic barriers. We shall discuss such periodic struc- 
tures with DEMVP barriers in detail in SccfVl 



IV. QUANTUM GOOS HANCHEN SHIFT IN 
SINGLE MVP AND EMVP BARRIERS 

In the preceding section, we have analyzed the trans- 
mission properties through a combination of electrostatic 
and magnetic vector potential barriers using the language 
of geometrical optics. In this section, we shall discuss 
another very important optical phenomenon, the Goos- 
Hanchen (GH) effec t 47 ' 48 , which is experienced by elec- 
trons incident on MVP and EMVP barriers. Recently, it 
has been shown that ballistic electrons passing through a 
p — n interface in graphene 4 ^ suffer a Goos-Hanchen shift, 
which changes sign at certain angles of incidence. Our 
treatment extends this analysis by including the effect of 
magnetic barriers on such a shift. 

The GH effect describes the shift in a beam of light 
suffering total internal reflection at an interface along 
the longitudinal direction (y-axis in our discussion). It 
has been known since the time of Newton^ 8 - and was first 
experimentally measured by Goos and Hanchen 4 ^. The 
shift is detectable since the extent of a real beam is al- 
ways finite. The shift occurs due to the fact that the 
totally reflected ray undergoes a phase shift as compared 
to the incident beam. The shift has an opposite sign 
if the second medium shows negative refraction and be- 
haves like a metamateria l 49 ' 50 . Such a lateral shift both 
for totally as well as partially reflected electron waves can 
also occur for non-relativistic electrons passing through 
a semiconductor barrier—, magneto-electric semiconduc- 
tor nanostructure^. In the subsequent discussion, we 
shall calculate the GH shift using the procedure given in 
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A(x) X t 




(b) 



<P" m* 



FIG. 8: GH shift for a wavepacket of massless Dirac fermions 
incident upon an EMVP barrier. The solid line corresponds 
to the upper component of the pseudospinor and the dotted 
line correspond to the lower component. The GH shift can be 
either (a) positive or (b) negative. 



terms by Taylor expanding around k and retaining only 
the first order terms. This yields 




qV/E p 



- -i 

I 



FIG. 7: The GH shift calculated for for single EMVP barrier 
of field strength (a) 0.1 Tesla and (b) 3 Tesla.The y-axis cor- 
responds to the incident angle The green region corresponds 
to the propagating wave solutions and thus the GH shift is 
set at 0. 



Ref4£ when massless Dirac fermions arc total internally 
reflected by MVP and EMVP barriers and analyze the 
results. We will also discuss how to calculate the GH 
shift for a DEMVP barrier and several such MVP and 
EMVP barriers as well. 

To calculate the GH shift, we consider the following 
wave packet (beam) of electrons impinging on an MVP 
or EMVP barrier exactly at the Dirac point: 



*m(x,y) 



ik y y-\-ik x (k y )x 



l 

e itf>(k v \ 



(37) 

The envelope function ensures the wavepacket is of fi- 
nite size along the y-direction and is sharply peaked at 
k y = k. Thus, k € (0,kp) and the angle of incidence 
(j>(k y ) € (0, f). Following standard procedure, we take 
the envelope Gaussian such that 



f(k y -k) = cxp[ 



(ky ~ ~kf 
2A^ 



(38) 



When Afc <C kp, we can approximate the k y dependent 



(j){ky) « (j){k) + —\- k {ky-k)] k X (ky) PS k X {k) + -^-\- k {ky- 

(39) 

Substituting the above results in Eq. ([3T)) and performing 
the Gaussian integrals, we obtain 



= j2TrA 2 k e l[ ' cv+k * Ck)x] 



e — T-IV-V+ r 
e -^lv-v'"] 2 e i<p(k) 



(40) 



where 



-k' x (k)x,yt 



-k' x (k)x-cl>'(k) (41) 



Thus, the upper and lower components of the spinorial 
wave function is localized at two different points along 
the j/-axis. 

The reflected wavepacket can also be written in an 
analogous way by making the transformation k x to — k x 
and 4> to 7r — 4> as well as multiplying the reflection am- 
plitude r(ky) = |r(fc y )|e i ^( fc «). The reflected wave can 
be straightforwardly written as 

/oo 
dk y f{ky - k)e lk 
-OO 

r(k y ) 



y — ik x (ky )x 



-i<p(k y 



(42) 



Here again s = sgn(E — V). For pure magnetic barrier 
this can be taken as 1. 

The spatial profile of the reflected wave function can be 
again obtained by first expanding all k y dependent quan- 
tities around k and retaining only the first order terms 
and then performing the Gaussian integrals in expression 
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(HI. This leads to 



with 



*r = 



^/2^2 e »[fcy-fex(S)x 
\r(k)\ 



e 2 
"-[y-V-] 



-7- i2 

v-v+] 



-i[<j>(k)-<t>' r (k)] 



(43) 



Here, y r + and y r _ are given by 



y\ = -<f>' r (k) + k' x (k)x, yL = -4>' r (k) + k' x (k)x + (/>'$) 

(44) 

Thus, the above expression shows that the upper as well 
as lower components get shifted because of the existence 
of the phase factor. The Goos Hanchen shifts of the 
upper and lower components are respectively given by 



y+ 
y r - 



y+ 
y™ 



-(j)' r {k) + 2k' x (k)x 
2<t>'(k)-ct>' r (k)+2k' x (k)x (45) 



Thus, the average GH shift for the electron falling on an 
MVP or an EMVP barrier can be taken as 



The situation is depicted with the help of a schematic 
diagram in Fig. [8] The average GH shift is not an origin 
dependent quantity. The last term in the above expres- 
sion is a coordinate dependent quantity and will get an 
equal and opposite contribution from the —<f)' r .(k) term 
upon the explicit evaluation of the reflection coefficient. 
The resultant expression will thus be independent of the 
choice of the coordinate of the interface from which total 
internal reflection will take place. Using the above ex- 
pression, one can calculate the GH shift when the angle 
of incidence <p is greater than the critical angle of inci- 
dence 4> c . 

For the case of either an EMVP barrier or an MVP 
barrier, the reflection coefficient r (k y ) can be calculated 
in the same way as done for the transmission coefficient 
given in Eq. (|16p by demanding the continuity of wave 
functions on both sides of the barrier interface at x = 
—d, and noting that on one side of the barrier the wave 
function is evanescent. Using the expressions used in Eqs. 
(HHJ) and (fTTj) , such a wave function can be written as 



a %k a x _|_ re ik^x x < _^ 

a , e - K{x+d) if x > _ d 



s U(k„x+<t,) _ re -i(k x x+<j>)} if x < _ d 



-i"/s , a , e- K ( x+d '> if x>-d 
Here, the constant 7 is given by 

K+ (ky — J-) 



(47) 



(48) 
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fCF 

K+ (ky — J^) 



MVP 
EMVP 



(49) 



k F , MVP 



(k y -— ) 2 -k 2 = k%, EMVP (50) 



where k' F has been defined in the Eq. (fP2]) . It may be 
mentioned that, for the wave packet considered above, 
these conditions for evanescent wave may not be satisfied 
for different k y that constitute the wave packet and only 
a fraction of such components may be totally reflected. 
However, as long as Afc <C k it is reasonable to assume 
that the entire wave packet is either partially transmitted 
or fully reflected. If, on the other hand, one considers a 
wave packet which is broader then the above conditions 
need to be relaxed. 

In view of the relations (fB"0"l) . it is instructive to 
parametrize k, k y 

in one of the following two alternative ways: 



-, for a totally reflected wave packet, 



\ky 



1 



— — — 1 = k' F (kp) cosha; k = k' F (kp) sinha, (E)MVP 



4>'(k)-cP' r (k) + 2k' x (k)x (46) or,fc^(fc F ) 



sin P; k 



cos/3, (E)MVP 



The second paramctrization can be heuristically under- 
stood as the angle made by the wavevector of the to- 
tally reflected wave with the surface normal in the first 
medium. With this notation, it can be verified that 
7 = exp(a) = cot(^). The continuity of the wave func- 
tion at x = — d gives the reflection coefficient as 



-ik x Dt 



J4> 



ss 7 



ss 7 



exp(—ik x D) exp(2i<5) (51) 



D that appears 
The contribu- 



Here tan<5 = tan</> + ss'j sec 4>, which is an expression 
similar to the one derived for scalar electrostatic barrier— 
with the exception of the prefactor e~ tk 
due to a different choice of the origin, 
tion to the GH shift from this extra term will cancel the 
contribution coming from the extra x dependent term in 
Eq. (|4"6"|) after setting x = —d. The reflection coefficient 
derived above is thus a unimodular complex number with 
the phase given by 

(j> r = -k x D + 25 (52) 
The expression of the GH shift will thus be 

cr = 4>'(k) - 25' (k) (53) 



where 



tan 1 (tan — A) 



Using this expression, one can calculate the GH shift for 
MVP and EMVP barriers when total internal reflection 
occurs. We shall next consider these cases separately 
with the help of representative figures and will also men- 
tion briefly how the GH shift can be obtained for a double 
EMVP or a periodic array of EMVP barriers. 
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1. GH shift for MVP barriers 



For pure MVP barrier, the total internal reflection will 
take place only when > <j> > — \ since the wave incident 
from the right and left hand side of the surface normal 
will behave differently^. Now for a given strength of the 
barrier the critical angle of total reflection is given by 

KplB 

Thus, there will be a finite GH shift for — <p c > (f> > — 5. 



2. GH shift for EMVP barriers 

For an EMVP barrier, total internal reflection can hap- 
pen for electrons incident from both sides of the surface 
normal but at different critical angles. For a given angle 
of incidence </>, it is possible to change the V adiabatically 
and get the electron wave reflected over a range of V sat- 
isfying | sin 6*| > 1. We can keep <j> fixed and increase V. 
At V = V c i the wave will be totally reflected such that 

sin \9\ = l,atV = V cl 



^ = l-[sin|0|-- 

hi kpiB 



],O<0< - 



= l-[ sm |0| + -J_],O>0>| (54) 

KplB 2 

Upon further increase of V the electron wave will re- 
main total internally reflected still the electrostatic po- 
tential reaches the second critical value of V C 2, such that 
sin \0\ = — 1. This gives 



,o<4>< 



1 + [sin |0| 



1 



W°^2 ^ 



For a given <fi for the range V c i < V < V C 2 the elec- 
tron wave will remain total internally reflected. This is 
demonstrated in Fig. [7] by plotting the GH shift over the 
entire range of V and <fr where the incident wave remains 
total internally reflected by the barrier. At all other re- 
gion in <j) — V plane the GH shift is set to ( green). At 
V = Ep, where the scale factor diverges, lies within this 
range. As discussed previously at this singular point the 
solutions of the wave equations correspond to zero mode 
solutions of the Dirac operator in presence of an EMVP 
barrier. These solutions in this case are always evanes- 
cent waves. Thus there will be finite GH shift at this 
point also which should be calculated after taking into 
account the special nature of the solutions here. Except 
at this singular point, at all other values of Using the ex- 
pression (|53[) the explicit expression for the GH shift for 



an EMVP barrier can be calculated from the following 
dimensionless expression 



1 - 2ff cos 2 5(1 + tan tan S + ss'^^] 



COS ( 



(56) 



The last quantity in the numerator can be equated 
to either coth a or sec (3 according to our parametriza- 
tion. As can be seen just at the point where TIR takes 
place this quantity is divergent. Accordingly the GH 
shift also diverges at this point. The calculated GH shift 
has been plotted in Fig. [7J for two different strength 
of the magnetic field. As one can see the border re- 
gion of finite GH shift, namely where the total inter- 
nal reflection takes place shows the highest value of GH 
shift. The lower and upper part of the curve again cor- 
responds to sin cf) — fc ^ < Oand > 0. Where as the left 
and right boundary corresponds to the cases of V < Ep 
and V > Ep. As can be seen in the lower part in most 
of the region the GH shift is negative with the left and 
right bordering region respectively showing the GH shift 
for positive and negative refraction. As explained the 
high value of GH shift just at the boundary where the 
total internal reflection begins either at a given <fi for a 
fixed V or at a given V for fixed <f> is due to the diverging 
behaviour of the term m the expression (f56|) . This 
can be related to the fact that just at the critical angle 
of the total internal reflection the wavevector lies in the 
interface of the two modeium. 

If we compare the above expression with the one cal- 
culated in Ref4& (Eq. 11), we shall see <j{4>) ^ — er(— </>) 
in our case. This happens because of the non-specular 
nature of electron refraction at an MVP— or EMVP bar- 
rier. This is an important difference in the quantum GH 
effect that occurs at the total internal reflection by MVP 
or EMVP barrier as compared to TIR by purely electro- 
static barrier. Such a quantum GH effect can lead to 
interesting devices in the regime of coherent electronic 
transport^. 

The GH shift at a double EMVP barrier the can be 
calculated in the same way. However here the wave suf- 
fers total internal reflection either at the first interface 
or at the second interface. Thus one needs to ascertain 
first whether the wave will be total internally reflected 
by the first or second regime. Once that is determined 
the calculation of the GH shift can be carried out in the 
same way as described in the case of single EMVP bar- 
rier. Once a wave suffers TIR it will be GH shifted and 
will not propagate further. However as we put an ar- 
ray of such barriers the fraction of the incident electrons 
suffering GH shift will be enhanced 



V. ELECTRONS IN A PERIODIC LATTICE OF 
DEMVP BARRIERS 

We shall now consider the motion of massless Dirac 
fermions in monolayer graphene through a periodic ar- 
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rangcmcnt of double EMVP barriers considered previ- 
ously in Section IIII Bl Recently, the transport of mass- 
less Dirac fermions through a periodic arrangement of 
pure magnetic barrier— and periodic as well as disor- 
dered superlattices of scalar potentials. To start with, 
we consider a finite number of repeated DEMVP barri- 
ers. In this case, the boundary conditions become im- 
portant. Such a structure forms a Bragg reflector as dis- 
cussed previously^. We then extend the analysis for an 
infinite lattice of such DEMVP barriers. 



A. Transmission through a finite number of double 
EMVP barriers 



structure with split gate voltage given in Eq. (fl7|) can be 
repeated infinitely. With a suitable voltage V, the split 
gate voltage can locally convert a charge-neutral region 
in monolayer graphene into ap — novn — p junction. 
Thus, such a periodic structure can also be thought of as 
a semiconductor heterostructure. As a result, the present 
problem will also address certain issues of magnetotrans- 
port through such a heterostructure when the magnetic 
field is highly inhomogencous and periodic. 

We consider each unit cell of size D = 2d for the MVP 
as well as for the electrostatic potential barriers. Thus, 
the n-th cell is given by (n — 1)D < x < nD. In the a-th 
part of the given unit cell, the wavefunction is given by 



Since the double EMVP barrier exhibits quite different 
transport properties as compared to a pure MVP barrier, 
it is interesting to consider how a Bragg reflector made 
of a sequence of double EMVP barriers will behave. This 
will also supplement our analysis about the reflection 
of such Bragg reflectors using double MVP barriers^!. 
The transmission through N such double EMVP barri- 
ers placed side by side can be calculated by generalizing 
the methodology described in the earlier section. The 
magnetic field for a Bragg reflector placed symmetrically 
around the origin is given byii 

B = B z (x)z = B£ B [S(x + nd) + S(x - nd) 

n-l 

+ (-l) p+n 2B5(x-pd))}z(57) 

p—l — n 

The series of wave function solutions in the various re- 
gions are linear combinations of right and left moving 
waves similar to the ones given earlier in Eqs. ([25]) and 
(f2"9")l for one double EMVP barrier. The transmission co- 
efficient can be calculated by using the continuity of the 
wavefunction in the interface of two such barriers. Just 
as Eq. (|3"5")) describes the solution for one double EMVP 
barrier, the solutions for n double EMVP barriers can be 
written in matrix form as 

[1 r] T = {M^) n {M s Mt)- r T» MVP M s ,Mt{M* A ) N [t of 

(58) 

Here Temvp is just the transfer matrix through a double 
EMVP barrier given in Eq ([55]). 

Representative plots are given in FigJSJa) for N = 5 
and in Fig^b) for N = 25 . The central column corre- 
sponds to the case when V = 0. As one increases |V|, the 
reflectance gets reduced and the transmittance increases. 
For low values of |V|, the reflectance is highly asymmet- 
ric due to the asymmetric behaviour of sin 9\ and sin 02 
with (j> (Fig(5j) . For higher values of |V|, full transmission 
occurs over the entire range of <fi. This significantly de- 
creases the reflectance of the Bragg reflector. However, 
with the increase in the number of double EMVP bar- 
riers this decay is slowed down and is preceded by an 
oscillatory behaviour of reflectance. 

B. Transmission through an infinite periodic lattice 
of DEMVP barriers 



h = a°e iq ^ x - nD) + by- iq ^ {x - nD) (59) 



d ■ s 1 

n ' n 



Here, a = 1,2; a* = a n ,b l n = b ni a 2 n = c ni b 2 n 
Si,s„ = S2',qi x = Qi,Q2x = 92- The exponential factor 
e ~nD revea j s the- existence of lattice translational sym- 
metry, which is not present for isolated single and double 
barrier structures discussed in Sections IIII Al and IIII Bl 
As compared to the case of pure MVP barriers, here q\ 
and q2 will change according to Eqs. (|2"l|) and (|2"B"|) . Also, 
Si 2 = sgn(_E ± V) can have same or opposite sign de- 
pending on the electrostatic potential V and the energy 
E. 

The continuity of the wavefunction at the first interface 
at x = (n — 1)D gives 



s 2 e 



1 1 

i&2 



S2e 

£ -iqiD 



-i8 2 



Cn-1 

d n -i 
b„e lqiD 



(60) 



Similarly, the continuity at the second interface at x 
(n — 1)D + d gives 



S2[e 



-i[q 2 S.-6 2] 



-s 2 e 



e «ji-- 

i[q 2 - 



(61) 



Using Eq. (j33"|) . Eqns. ([60|) and (fBTj) can be rewritten as 



M S2n _ 1 M 02 



Cn-l 
d n -l 



M Sln M 9l M Bl 



M slin M Sl M Bl 

M s Mg Mb 2 



(62) 



The above two matrix equations can be combined as 



In this case, to simplify the analysis, we shall ignore the 
boundary effects and we assume the double MVP barrier 
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B=3T, qV/^-=-4.2 B=3T, qV/E^-1 .4 B.3T, qV/Ep.O B.3T, qV/E^I .4 B=3T, qV/^-=4.2 




B=3T, qV/Ep=-4.2 B.3T, qV/^-~ 1 .4 B.3T, qV/E^O B=3T, qV/E^I .4 B=3T, qV/Ep-4.2 




FIG. 9: Transmittance and Reflectance through a finite sequence of 5(upper three rows) and 25 (lower three rows) DEMVP 
barriers. On the x-axis is plotted the incident angle <j>. 



Cn-l 
d n -i 



Mq 2 1 M~ 1 M S1 Mg 1 M Bl 1 M~ 1 M s2 Mq 2 Mb 2 

I 



However, according to Bloch theorem, 



C n -1 


= e- iKD 











(64) 



The matrix K„ 



is thus unimodular. 



if li 

K 2 i K22 

Combining the results in Eqs. (j6"3")) and (f6~4"|) we obtain 
the eigenvalue equation 





K 12 ~ 






= e~ iKD 


Cn-l 


K21 


K 2 2_ 




d„_ 




d n -l 



where K is the Bloch momentum. The explicit form of 



c» 




'K u 


Ki2~ 




c„ 


d„ 




K21 


K 2 2_ 




d n _ 



(63) 



the various terms in the matrix can be written as 
#ii = F + {8 2 ,8 1 )F+{8 ll 8 2 )e-' l ^ + ^ d 

+F-(e 2 ,e 1 )F-*{e 1 ,e 2 )e i ^-^ d 

K 12 = F+(d2,8 1 )F-(8 1 ,8 2 )e-^-^ d 
+F-(6 2 , e 1 )F + *{e u 92)e^ qi ~ q ^ d 

K21 = F-*(62,8i)F+(6 u 62)e- l ^ + ^ d 
+F+*(82,8 1 )F-*(9 1 ,02)e t ^-^ d 

K22 = F + *(8 2: 8 1 )F*- (0 1 ,0 2 )e i ^-^ d 

+F-*(8 2 , 6 1 )F+(6 1 ,62)e- l{qi+q ^ d (66) 



where 



F ± (6 k ,9i)=e- i6k ±s 1 s 2 e iei ,toT k,l= 1,2. 

(65) 

The eigenvalues A which are complex conjugate are given 
by 

det\K mat -XI\ = 

^Ai+A 2 = exp(-iKD)+exp(iKD), (67) 
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B=3T, qV/E F =-4.2 

A 



If ^ 



B=1T, qV/E F =-4.2 
i 




B^O.IT. qV/^_=-4.2 



1 



B=3T, qV/E^-1.4 




B=1T, qV/E^-1.4 



V _ 



B=0.1T, qV/E^-1.4 



B=3T, qV/Ep-0 



B=1T, qV/Ep-0 




B=3T, qV/Ep-1 .4 



B=1T, qV/Ep-1.4 



B=0.1T. qV/E^I.4 



B=3T, qV/Ep-4.2 



u 

B=1T, qV/^4.2 




-71/2 -71/4 71/4 7I/2-7I/2 -71/4 71/4 7I/2-7I/2 -7D/4 71/4 7C/2-7I/2 -71/4 71/4 7t/2-7l/2 -71/4 71/4 71/2 



FIG. 10: The allowed and forbidden transmission regions for EMVP barriers of strengths 0.1, 1 and 3 Tesla. 



which finally gives 



K(d>,B) 



2d 



cos- 1 [-Tr(^)] 



(68) 



The condition \iTr(Kij)\ < 1 corresponds to propagat- 
ing Bloch waves whereas \^Tr(Kij)\ > 1 leads to evanes- 
cent Bloch waves that correspond to forbidden zones 
in the band structure in presence of periodic MVP + 
electrostatic barriers. Writing explicitly in terms of the 
wavevectors q\ , q2 and the angles 01,02, the above eigen- 
value condition reads 



cos KD 



cos q\d cos q2<i + sin q\d sin q2(L x 
[ tan 0\ tan 62 ' 



r 2 „ ] (69) 

cos 0i cos 02 



The above equation provides the band structure for 
a periodic EMVP barrier. This is an extension of 
the Kronig-Penney (KP) model to the case of two- 
dimensional masslcss Dirac fermions. Thus, it is inter- 
esting to compare the band structure that arises out 
of this model and discuss its differences with the band 
structures of other variants of the KP model. The origi- 
nal KP model describes the transmission of Bloch waves 
through a one dimensional periodic potential^. Several 
authors have also studied the relativistic version of the 
KP mode l 55 ' 56 i 57 ' 58 . In these studies, the motion con- 
sidered is strictly one dimensional. The non-rclativistic 
KP model in periodic structures created by MVP barri- 
ers has also been studie d 37 i 39 . This second set of studies, 
being for non-relativistic electrons, has a different set of 
boundary conditions at the unit cell interfaces since the 
generic wave equation in that case is a second order equa- 
tion. The current problem is a relativistic KP model for 
two dimensional masslcss fermions in the presence of a 
potential that is effectively in one dimension and is then 



expected to show souk- differences as compared to the 
previous ones. For electrons in graphene, such problems 
have been studied very recently for periodic structure of 
different types of magneti o 17 i 18 as well as electrostatic^ 
barriers. The novelty in the present problem is that it in- 
cludes the combined effect of electrostatic and magnetic 
vector potentials. In the following analysis, we will see 
that it adds some interesting features to the band struc- 
ture. 

In FigfTU] is plotted the quantity cos K D given by 
Eq. (|6"9"|) as a function of <j> for three different values of 
the magnetic field. The central column corresponds to 
the case of V = 0, the case of an infinite series of pure 
magnetic barriers. The left and right columns plot the 
same quantities for a range of voltages with ± sign. At 

V = 0, the plot is symmetric about = 0, which also 
implies the symmetric transmission property of a dou- 
ble MVP barrier structure^ 7 , which forms the unit cell of 
such a periodic structure. Also, a forbidden region opens 
around = and widens with increasing barrier strength 
B. 

As we change the value of V = ±0.1V (equivalent 
to zkqV/Ep = 1.4), the allowed region of propagation 
shrinks as compared to the case for V = at finite val- 
ues of |V|. In addition, the distribution of the allowed 
and forbidden regions becomes asymmetric on the 4> axis. 
Beyond a certain angle, which is different for positive and 
negative segment of the 0-axis, for V = 0.1V as well as 

V = — 0.1V there is no transmission since the entire re- 
gion is forbidden. When V is increased further in mag- 
nitude to near V = IV (equivalent to ztqV/Ep = 4.2), 
for one sign of </> all solutions correspond to propagat- 
ing Bloch waves whereas for the segment of 0-axis that 
corresponds to opposite sign of <p such propagating so- 
lutions exist only over a part. Upon changing the sign 
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of the V, the behaviour in these two parts of the 0-axis 
get exchanged. This asymmetric behaviour as a function 
of cf), as well as the reversal upon the change of sign in 
V, persist for other values of B. However, when B is in- 
creased from B = 0.1T to B = 3T, more forbidden zones 
form in between the conducting regions on the 4> axis. 
All these observations are consistent with the plots of 
refraction angle and transmission for the double EMVP 
barrier given in Fig|5]and Fig(6j 



B-3T 




qVIE F 



FIG. 11: Band structure in presence of periodic DEMVP bar- 
riers: The colour axis plots the energy in the (k y , V) plane 
for barriers of strengths (a)0.1 Tesla and (b)3 Tesla. White 
regions are the forbidden zones. 

In Figfrrjis plotted the full dispersion relation by vary- 
ing energy E. Since the discussion treats the quasipar- 
ticles as massless Dirac fermions, our treatment is valid 
as long as the energy E is close enough to the Dirac 
point such that the linear dispersion relation between 
the energy and momentum is maintained. Accordingly, 
in FigJTSJ a portion of FigfTT] for certain characteris- 
tic values of V has been magnified and presented. In 
this regime, for a given value of energy, we can write 
E = hvp\k\ and k y = |fc|sin<^. Thus, for a given E , 
\k y \ G {7^,0} and \<j>\ G {§ ,0}. Given the values of 
cj), it is possible to determine with the help of Eq. (j6"51) 
whether the Bloch waves will be propagating or evanes- 
cent, the latter case corresponding to the forbidden zone. 
Using the above criterion, we can plot E in the plane of 
k y and V and identify the forbidden regions. This gives 
us the band structure. The wavevector k y is plotted units 
of j- . The corresponding plots has been displayed in two 
figures over a range of V for two different values of mag- 
netic field in Fig. 1111 The lower plot corresponds to a 



field strength of B = 0.1T whereas the upper plot corre- 
sponds to the field strength of B = 3T. Since the unit of 
k y , that is (VTb) -1 is higher in the case of B = 3T, the 
upper limit of y axis is different in these two plots. We 
shall mostly concentrate on analyzing the lower part of 
the diagram where the Dirac fermion description is valid. 

The plots show conducting regions over a wide range of 
energy with forbidden regions in between. For B = 0.1T, 
the region near V — is mostly conducting for different 
values of energy where the forbidden region starts open- 
ing up in both left and right directions of the V axis. 
Let us emphasize again that exactly at = 1, the 
above method of obtaining the band structure through 
Eq. ([69]) is not valid since the solutions of the wave 
equations are always evanescent waves as explained in 
section UlI Al However we can try to analyze the behav- 
ior close to this point within this method. Near j^. — > 1, 
an extended forbidden zone appears. To understand the 
behaviour in this region more explicitly we have plotted 
in Fig[T2]the band structure at = 0.7V which is close 

to the singular region — > 1. For this value of V, the 
conducting regions are intercepted by large patches of 
forbidden zones. However, if the voltage is switched to 
a negative value, namely at ^ = —0.7, the behaviour 
is completely changed. Now the forbidden region in the 
band structure over the same range of kyig shrinks con- 
siderably (left column of the figure) higher transmission. 
This behaviour can be directly related to the laws of re- 
fraction for a double EMVP barrier given by Eqs. (pH})) 
and i[3l]). 

In Fig. rn] for higher values of V where > 1, the 
system is conducting for allmost all values of E. This is 
expected also from the earlier analysis (FigO and FiglB]) 
of such transport through one DEMVP barrier. For a 
higher magnetic field of B = 3T, near V = a forbidden 
region opens up at various values of the energy. Also the 
forbidden region near |^ — * 1 is much larger as com- 
pared to those for B = 0.1T. This can be seen from 
Fig. [TT]as well as FigfTJ] These figures also suggest that 
the gapped regions to the left and right of V = are lo- 
cated in a pronounced asymmetric manner as compared 
to the case when B = 0.1T. We may recall that such 
pronounced asymmetry in the transmittancc at higher 
values of B = 3T was also seen in the plot of transmit- 
tance through one single double EMVP barrier in Fig(S] 
This asymmetry in the band structure of transport as a 
function of the applied voltage as well as the opening up 
of large forbidden zones for certain values of V provides 
a way to tune the transport properties of massless Dirac 
fermions in monolayer graphene in a controllable manner. 

VI. CONCLUSION 

In the above analysis, we have considered in detail 
the transport of massless Dirac fermions in graphene 
through a combination of singular magnetic barriers and 
applied electrostatic voltage. Within a coherent ballis- 
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(a) B = OAT, qV/E F =-0.7 
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FIG. 12: Band structure at some characteristic values of V with B-field strength strengths (a) 3 Tesla, and (b) 0.1 Tesla. 



tic transport regime, the problem also addresses the ef- 
fect of highly inhomogeneous magnetic field on transport 
through a graphenc-bascd p-norn-p type junction as 
well as heterostructures consisting of such junctions. We 
have shown how the transmission through such barriers 
gives rise to several analogues of interesting optical phe- 
nomena such as total internal reflection for positive as 
well as negative refraction and a quantum version of the 
Goos Hanchcn shift. We have shown that with the help 
of this optical analogy established for single or double 
EMVP barriers, all the major characteristics of the band 
structure in the presence of a series of such EMVP bar- 
riers can be well understood. All these calculations have 
been done in the close vicinity of Dirac point and with 
the assumption that the K and K' valleys are degener- 
ate. Thus, the conclusions are applicable in the regime of 
coherent ballistic transport. Present experimental tech- 
niques render such a regime accessible and their scope 
is widening further—. To explore the full band struc- 
ture, one needs to calculate the impact of such periodic 
EMVP barriers within the tight binding approximation. 
However, while doing such tight binding calculations, it is 
important to note that the lattice spacing in graphene is 
generally much smaller as compared to the typical width 
of such highly localized magnetic barriers (of the order 
of a few nanometers). Such calculations must take into 



account the finiteness of the barriers. A comparison with 
other works carried out for periodic structures with finite 
magnetic barriers^ shows that many properties related 
to the transmission through such barriers are very similar 
to the ones we have obtained by treating this barrier in 
the delta function approximation. Finite size effect and 
boundary conditions can strongly affect the conclusions 
based on the above treatment^. It will be interesting to 
consider finite size effects on the above mentioned band 
structure properties by carrying out similar calculations 
for various types of graphene nanoribbons. To summa- 
rize, if one is interested in the low energy transport prop- 
erties around the Dirac point, the preceding calculations 
provide a number of interesting results which we believe 
to be useful for graphene based electronic s 60 i 61 . 
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